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Commonly Observed Molecular Dense Gas Tracers 
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ABSTRACT 

The optically thin critical densities and the effective excitation densities to 
produce a 1 K krn/s (or 0.818 Jy km/s ( 10 qq H z ) 2 (%7f 1 ) 2 ) spectral line are tab¬ 
ulated for 12 commonly observed dense gas molecular tracers. The dependence 
of the critical density and effective excitation density on physical assumptions 
(i.e. gas kinetic temperature and molecular column density) is analyzed. Critical 
densities for commonly observed dense gas transitions in molecular clouds (i.e. 
HCN 1 — 0, HCO + 1 — 0, N 2 H + 1 — 0) are typically 1 — 2 orders of magnitude 
larger than effective excitation densities because the standard definitions of criti¬ 
cal density do not account for radiative trapping and 1 K km/s lines are typically 
produced when radiative rates out of the upper energy level of the transition are 
faster than collisional depopulation. The use of effective excitation density has a 
distinct advantage over the use of critical density in characterizing the differences 
in density traced by species such as NH 3 , HCO + , N 2 H + , and HCN as well as their 
isotpologues; but, the effective excitation density has the disadvantage that it is 
undefined for transitions when E u /k T for low molecular column densities, 
and for heavy molecules with complex spectra (i.e. CH 3 CHO). 


1. Introduction 

The density at which a particular molecular transition is excited has always been an 
interesting quantity for studies of molecular gas in the interstellar medium. Two concepts 
have been developed, the critical density and the effective excitation density, to quantify the 
excitation density of a transition. A traditional approach in studying the critical density has 
been to consider the molecular rotational energy levels as a simple two level system. While 
this approach has pedagogical advantages, the interpretation of concepts, such as critical 
density, are oversimplified and incorrect. Molecules are fundamentally multi-level systems 
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with radiative and collisional processes between the levels (Figure 1). It is almost never 
appropriate to approximate commonly observed dense gas tracers, such as HCO + , HCN, 
NH 3 , N 2 H + , etc., as two level systems. Therefore, in this tutorial, we eschew the two level 
approximation and consider the full multi-level nature of molecular excitation. It is assumed 
that the reader has a basic familiarity with statistical equilibrium, radiative transfer, and 
molecular spectroscopy. For a summary of radiative transfer techniques and molecular data 
used in those calculations, see van der Tak (2011) and references therein. The purpose of 
this tutorial is to quantify and systematize the use of critical density and effective excitation 
density for dense gas molecular tracers. 

We shall use the excitation of two molecules, HCO + and NH 3 , as examples throughout 
this tutorial. Figure 1 shows the energy level structure for both molecules. HCO + is an 
abundant linear molecule with a simple electronic ground state meaning that there is 
zero net electronic spin and zero net electronic angular momentum that could interact with 
and split the rotational levels of the molecule. Electric dipole selection rules limit transitions 
to A J = ±1. The 1 — 0 ground state rotational transition is a millimeter transition at 89 GHz 
(3.3 mm) and the frequency of rotational transitions progress linearly into the submillimeter 
part of the spectrum. 

NH 3 is an abundant pyramidal symmetric top molecule with a simple 1 A\ electronic 
ground state also meaning that there is zero net electronic spin and zero net electronic angular 
momentum. The energy levels are designed by Jk where J is the total rotational angular 
momentum and K is the projection of the angular momentum onto the NH 3 symmetry 
axis. Electric dipole transition do not allow the projected angular momentum to change 
(A K = 0). As a results, the energy level structure is usually organized in K-ladders with 
each K-ladder designated as either ortho (K = 0, 3, 6, etc.) or para (K = 1, 2, 4, 5, 
etc.) due to the spin symmetry required to satisfy Fermi-Dirac statistics for each rotational 
energy level. The plane of three H atoms can tunnel from positions above and below the 
N atom resulting in inversion splitting of each rotational level (with symmetry denoted by 
±) except for levels with K = 0 for which half of the inversion levels are missing due to the 
Pauli exclusion principle (see Townes & Schwalow 1975 for a detailed explanation). Electric 
dipole transitions are allowed between the inversion-split levels with different symmetries and 
these inversion transitions occur at wavelengths near 1.3 cm (24 GHz). Since J and K do 
not change in a pure inversion transition, the shorthand notation (J, K ) has been adopted. 
As with HCO + , there are rotational transition with A J = ±1, but since NH 3 is a lighter 
molecule than HCO + ( E rot oc 1/J where I is the moment of inertia of the molecule measured 
from a coordinate system at the center-of-mass), the transitions occur at submillimeter and 
far-infrared wavelengths. Thus, the lowest energy transitions of HCO + and NH 3 cover a 
wide range of wavelengths from the centimeter to the far-infrared. 
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Fig. 1.— Left: The energy level structure for HC0 + showing the lowest 4 energy levels. 
Electric dipole allowed transitions are shown along with the quantum numbers (J) to the 
right of each level and statistical weights in red to the left of each energy level. Right: The 
energy level structure for Nff 3 showing the lowest 9 energy levels. Electric dipole allowed 
transitions are shown along with the quantum numbers (Jr) to the right of each level and 
statistical weights in red to the left of each energy level. This figure was adapted from figures 
in Mangum & Shirley (2015). 
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2. Critical Density 


The critical density has traditionally been used as a measure of the density at which a 
particular transition is excited and is observed at radio wavelengths. The definition of the 
critical density is not consistent throughout the literature. Some definitions only consider the 
two energy levels involved in the transition (2 level approximation) while other definitions 
use the multi-level nature of collisions to sum over all collisions out of the upper energy level 
or only from the upper energy level to lower energy levels. All of the standard definitions 
assume the molecular emission is optically thin and ignore radiative trapping which is very 
important and cannot be ignored for many commonly observed dense gas tracers (e.g. HCN 
1 — 0, HCO + 1 — 0, N 2 H + 1 — 0, etc.; see §2.2). The standard definitions also usually ignore 
continuum backgrounds (see §4). 


We start by defining the critical density as the density for which the net radiative decay 
rate from j —> k equals the rate of collisional depopulation out of the upper level j for a 
multilevel system. The net radiative decay rate from j —> k is the spontaneous decay rate 
from j —y k plus the stimulated emission rate from j —>■ k minus the absorption rate from 
k j. Mathematically, the ratio of rates is 

_ n jAjk + n jBjk'U>v jk ~ '^kBkjUy jk 
n c n j Iji 

where n 3 is the number density (cm -3 ) of molecules in the upper level of the transition, is 
the number density (cm -3 ) of molecules in the lower level of the transition, n c is the number 
density (cm -3 ) of colliding particles (typically H 2 , H, or e - ), and j 3l are the collision rates 
(cm 3 s -1 ) out of level j into another level i. n c is equal to n crit when 9ft = 1. Ajk and B 3 k 
are the Einstein A (s -1 ) and B (erg -1 cm 3 s -2 ) coefficients related by 


9k R _ R _ c 

— a kj — — 

9j 


8 vhtjA* 


( 2 ) 


u Vjk is the energy density (erg cm -3 Hz -1 ) of the radiation fields at the frequency of the 
transition Ujk. The energy density is written in terms of a Planck function and is related to 
the angle and polarization averaged photon occupation number, n p h, by 


u v {T) 


87 Thu 3 1 

£.3 g hv/kT _ ^ 


87r/iz/ 3 


n ph (T , v) 


( 3 ) 


The total energy density can include contributions from the cosmic microwave background 
(CMB), dust continuum emission from the interstellar radiation field, dust continuum emis¬ 
sion local to the molecular cloud, optically thick molecular emission at the frequency z/^, 
and other continuum and spectral line sources. 
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2.1. Optically Thin Approximation 


We define the optically thin critical density without a background by solving for n c in 
Equation 1 when 3ft = 1 and ignoring the stimulated emission and absorption terms: 


n 


thin,no bg _ 

crit 


A 


jk 


Ajk 


%/j fO' i y - / j<j t y 'j 


i>j 9j 


9±ry. ,p—(Ei—Ej)/kT k 

hjc. 


( 4 ) 


The sum over collision rates in the denominator of Equation 4 has been split in two terms 
depending on whether the collision from level j is to an energy level lower in energy (first 
sum) or higher in energy (second sum) than level j. Upward collision rates are related to 
downward collision rates using detailed balance by 


- ry..9i Le -{E i -E j )/kT k 

jji hj c 

9j 


for Ei > Ej 


( 5 ) 


where T k is the kinetic temperature of the collisional partner. Published collision rates (i.e. 
the Leiden molecular databas^; Schoier et ah 2005) quote only downward collision rates. 
In the two level approximation, which is widely quoted but rarely appropriate, the sums 

thin,no bg Aj & a i , thin,no b 

crit'2lvl ‘ — ^7’ AS a resui U n crit,2lvl 


i i 1 j thin,no bq 

reduce to a single term such that n 


■- —. As a result, always be larger 

'Yj fc 

(usually by at least a factor of a few) than n^’ n ° b9 ■ 


It is instructive to compare how the optically thin critical density compares to the exci¬ 
tation temperature. For a particular transition between levels j —>■ k, the level population are 
related to the excitation temperature, T ex , between the two levels by Boltzmann’s equation 


Tij gj 

— = — exp 

W'k 9 k 


Ej — E k 

kT P r 


for Ej > E k 


( 6 ) 


The excitation temperature (and therefore the level populations) depend on the density of 
colliding partners (n c ), the gas kinetic temperature of colliding partners (T k ), and the energy 
density of radiation fields that the molecules are exposed to. The dependence of T ex on the 
density and gas kinetic temperature for the ground state transitions of ortho and para NH 3 
are shown in Figure 2. The only background radiation field assumed in this example is 
the CMB. At low densities (n c < 10 3 cm -3 ), the level populations come into equilibrium 
with the CMB and T ex approaches T cm b = 2.725 K. As the density increases, collisions 
become more important in determining the level populations. The excitation temperature of 
(1,1) centimeter wavelength inversion transition equilibrates to T k by densities > 10 5 " 5 cm -3 
while the submillimeter 1 0 — 0q rotational transition equilibrates to T k at significantly higher 
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Fig. 2.— Top: The excitation temperature is plotted for the NH 3 (1,1) inversion transition 
(23.7 GHz, 1.3 cm) and 1 0 — Oo rotational transition (572.5 GHz, 0.52 mm) for a gas kinetic 
temperature of 10K and a total NH 3 reference column density of logiV = 14.8 cm -3 (an 
ortho to para ratio of 2 : 1 is assumed). The blue dashed line and the green dotted line 
indicate n^’ nob9 for the (1,1) and 1 0 — 0 0 transitions respectively. Bottom: The variation of 
excitation temperature for the NH 3 (1,1) transition with gas kinetic temperature is plotted 
for Tk = 10, 20, and 30 K. 
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densities (> 10' cm -3 ). The general shape of the T ex with logn c curves is a sigmoid curve 
with stimulated coupling with the background radiation field setting the lower bound at low 
densities and collisions with colliding partners at T k setting the upper bound. At densities 
where T ex < T k , the level populations are said to be sub-thermally populated. This is the 
situation for most of the commonly observed dense gas transitions in the ISM. When T ex 
approaches T kl the level populations are said to be thermalizedE 


At optical wavelengths (the high frequency limit), critical density is interpreted as the 
density at which an atomic forbidden line is quenched by collisions meaning that a collisional 
de-excitation occurs before a photon can be generated from spontaneous emission. In the 
optical, quenching occurs when T ex is thermalized ( T ex —> T k ). At radio wavelengths (the low 
frequency limit) the critical density occurs at different positions along the sigmoid T ex curve 
for different molecular tracers and different transitions. In the case of the 1.3 cm (23.7 GHz) 
NH 3 ( 1 , 1 ) transition, the optically thin critical density is two orders of magnitude below the 
thermalization density (Figure 2). Early radio astronomy observations of molecules such as 
OH and NH 3 occurred at centimeter wavelengths and since the critical density at these low 
frequencies is just above the regime where the T ex starts to rise above equilibration with the 
CMB, critical density was used to indicate the density at which a transition would appear 
in emission. As the frequency of the transition increases, the critical density moves farther 
up the T ex curve. For the 0.52 mm (572 GHz) NH 3 lo — Oo transition, the critical density 
is nearly equal to the density required for thermalization. Millimeter transitions are an 
intermediate case; thus, neither the radio interpretation of critical density being the density 
at which an emission line is excited nor the optical interpretation of critical density being 
the density at which an emission line is quenched is appropriate. 


2.2. Optically Thick Emission 

In reality, radiative trapping is very important for many molecular transitions observed 
in the dense interstellar medium and cannot be ignored. The energy density is then modified 
using the solution to the equation of radiative transfer to include emission from the molecule 
at the excitation temperature describing the level populations for a transition from j —» k 


2 Note that the T ex curve for NH 3 (1,1) shows a region of super-thermal excitation ( T ex > T k ) at densities 
immediately before thermalization for the T k = 30 K curve. This is due to a bottleneck in radiative 
transitions, caused by the five orders of magnitude difference between the spontaneous emission rates for 
the inversion transition and the THz rotational transition into the upper inversion level, resulting in a slight 
overpopulation of the upper level of the inversion transition. This effect is more pronounced for higher T k 
because upward collisions to the Jk = 2i levels are more efficient at high T k (see Equation 5). 



as 


Uv Jk = u ^k( T bg)ft + u Ujk (T ex ){ 1 - ft) , (7) 

where ft is the solid angle averaged escape fraction (see Chapter 19 of Draine 2011 for a 
derivation). If we multiply the energy density by the conversion factor between Einstein Bjk 
and Ajk terms from Equation 2, then we find that 

c 3 - - 

87r/lZ/ 3 blphiThg, Vjk)P T blph(Texi b'j /c) (1 ft') ■ (8) 


In the limit that the background radiation can be ignored ( u Vjk {T bg )ft <C u Vjk (T ex )(l — ft)), 
the critical density becomes 


thick,no bg 
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n 3 T?'i 


Tlj Aj 
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n 3 Tji 

1 + ( n p h(T e l,u jk )) n Ph( T ex, Vjk)(ft ~ ft) 


E*,- T? 


ft Aj k 


( 9 ) 


Comparing Equation 9 with Equation 4, we see that the effect of radiative trapping is to 
reduce the spontaneous transition rate by the escape fraction ft Ajk (see Scoville & Solomon 
1974, Goldreich & Kwan 1974). 

The relationship between ft and the line of sight optical depth, r Ujk , depends on the 
geometry and kinematics of the region. In the case of a static, uniform density sphere 


ft = 


4r ^ 8r 4 


+ e 2T ^k 


3 3 

+ 


4t 2 8 r 3 


( 10 ) 


(Osterbrock 1989), while in the case of a spherical cloud with a large velocity gradient v oc r 
(the LVG or Sobolev approximation), 


ft = 


(1 - e~ Tv i k ) 


r ^fc 


in) 


(Sobolev 1960, Castor 1970). For large optical depths {r Ujk >~ 5), the effective spontaneous 
transition rate is then given by ft Ajk ~ Ajk/r Ujk . 
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This result has an important effect on the critical density by reducing its value from 
the optically thin value. The most commonly observed transition of molecular gas is 12 CO 
1 — 0. The optically thin critical density at T k = 10 K is ri^’ no bfJ = 1 x 10 3 cm -3 ; however, 
since 12 CO 1 — 0 has optical depths that are typically more than a factor of 10 in molecular 
clouds (and usually much higher), it is easy to observe strong lines (> 1 K krn/s) in gas that 
has densities < 10 2 cm -3 . Optical trapping is so important for CO that it does not make 
sense to use the optically thin critical density. 

Commonly observed lines of dense gas tracers in cores and clumps also tend to be 
optically thick. The upper left panel of Figure 3 shows the predicted optical depth of HCO + 
1 — 0 through 4 — 3 at different densities from a simple (single density, single temperature) 
radiative transfer model predicting that these transition are expected to be very optically 
thick over a wide range of densities. Observations support this prediction. For example, 
several hundred high-mass clumps mapped in the MALT90 Galactic survey have a median 
optical depth of r = 23 for the HCO + 1 — 0 transition (Hoq et al. 2013). Clumps observed 
in the Bolocam Galactic Plane Survey (BGPS) in the higher frequency transition HCO + 
3 — 2 have an average r = 10 (Shirley et al. 2013). As a result, the analytically calculated 
critical densities using the optically thin approximation (Equation 4) are systematically 
biased toward higher densities by an order of magnitude or more. HCO + is not unique 
as most commonly observed dense gas tracers (i.e. HCN, CS, N 2 ID, etc.) in the dense, 
cold neutral medium are typically both sub-thermally populated ( T ex < T k ) and optically 
thick (also see Plume et al. 1997, Gerner et al. 2014). In practice, to properly estimate 
the optical depth, a coupled radiative transfer and statistical equilibrium calculation must 
be performed to find the level populations upon which the optical depth depends (since 
r Ujk oc f(n k B k j — rijBj k ) ds along the line-of-sight s). Since the escape fraction depends 
nonlinearly with column density (for modest optical depths), nonlinearly with geometry, 
and nonlinearly with gas kinetic temperature, it is difficult to make a simple tabulation of a 
p correction to n crit y . 


3. Effective Excitation Density 

It is possible to detect strong (1 K krn/s) molecular lines at densities well below the 
critical density (see Figure 3). An alternative definition has been developed to trace the 
effective density at which a modest line intensity is observed that is based on radiative 
transfer calculations with reasonable assumptions about the column density and gas kinetic 
temperature of the region (Evans 1999). The effective excitation density is defined as 


^e// n c 3 


T r dv = lKkm/s , 


( 12 ) 
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Fig. 3.— Top Left: Line of sight optical depth for HCO + transitions. A column density 
of log IV = 14.0 and a gas kinetic temperature of T k = 10 K were assumed. Top Right: 
The integrated intensity for ffCO + transitions. A column density of log N = 14.0 and a gas 
kinetic temperature of T k = 10 K were assumed. The horizontal dashed line is the criterion 
for the definition of n e ff used in this paper. Bottom Left: n e // is inversely proportional 
column density. Only the ffCO + 1 — 0 and 4 — 3 transitions are shown for clarity. The 
solid line is calculated at T k = 20 K while the dashed lines above and below the solid line 
correspond to T k = 10 K and T k — 30 K. Bottom Right: The dependence of n e ff on choice 
of gas kinetic temperature for the four lowest rotational transitions of ffCO + . n e ff was 
calculated assuming the reference column density of log N = 14.0. 
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the density which results in a molecular line with an integrated intensity of 1 K krn/s. The 
choice of 1 K krn/s is arbitrary but represents an easily detectable integrated intensity IE For 
lower integrated intensity criteria, n e ff decreases (Figure 3). It is important to note that the 
effective excitation density is based on a simple empirical criterion and is not based on 5ft = 1 
or any specific ratio of radiative rates to collisional depopulation rates. n e // does however 
account for the effects of radiative trapping. The original definition presented in Evans (1999) 
defined n e ff as the density resulting in T R — 1 K line and used a LVG radiative transfer 
code to calculate n e /f. We shall use the radiative transfer code RADEX (van der Tak et al. 
2007) with the escape probability formalism for a static, uniform density sphere, with the 
CMB as the only background radiation held, and with H 2 as the dominant collision partner 
for all calculations. The definition in Equation 12 is approximately equivalent to the Evans 
1999 definition when a FWHM linewidth of 1 krn/s is assumed in RADEX calculations. 


3.1. Properties of Effective Excitation Density 

The column density of the observed molecule and the gas kinetic temperature of the 
colliding partner must be assumed in order to calculate n e //. The sensitivity of n e // to these 
assumptions is explored in Figure 3 for the HCO + molecule. At lower column densities, a 
higher n e // is required to produce a 1 K krn/s line. As the column density increases by 
an order of magnitude, the n e ff decreases by almost exactly an order of magnitude. The 
effective excitation density is inversely proportional to column density 

N f 

n eff(Tk, N ) ~ n e ff(Tk, N re f ) , (13) 

where N re f is a reference column density. This simple scaling proportionality breaks down 
for n e ff >5 x 10 6 cm -3 and the relationship between n e // and N becomes noticeably non¬ 
linear (see the T k = 10 K curve for HCO + 4 — 3 in the bottom left panel of Figure 3). In 
those cases, a new RADEX model should be run to scale tabulated results to a new column 
density. 

The dependence of n e ff to T k is nonlinear and more sensitive at low kinetic temperatures 
(Tfc < 20 K) due to the inefficiency of upward collisions at lower T k (Figure 3). An order of 
magnitude higher n e ff is required to make a 1 K km/s line for T k = 10 K than at T k = 30 
K for HCO + 3 — 2. This sensitivity is due to the exponential dependence from upward 


O . 2 

J The radiation temperature is related to the observed flux density S v , k through Tr = R 0 - S v . k 

^‘ Kl/ jk^ l ‘^ earn J 

n0 2 

where Qbeam = 41 *^ is the solid angle of the telescope beam with FWHM Qbeam • AIK km/s integrated 
intensity corresponds to an integrated flux density of 0.818 Jy km/s ( iooghz ) 2 ) 2 - 
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collisional rates (Equation 4) and becomes more pronounced when the upper energy level 
is significantly higher than T\ (i.e. the effect is stronger for 4 — 3 vs. 1 — 0 at T k < 20 K 
because E u /k = 42.8 K for J = 4 vs. E u /k = 4.3 K for J = 1). 

Table 1 gives n e // for 12 dense gas tracers at T k = 10, 15, 20, 50, and 100 K for a reference 
column density. Recent surveys of dense gas in the Milky Way may be used to guide our 
choice of the reference column density. For example, observations of several hundred clumps 
by the MALT90 (Meittenen et al. 2013, Hoq et al. 2013), ChaMP (Barnes et al. 2012), 
and BGPS (Shirley et al. 2013) surveys indicate that HCO + column densities typically span 
from logiV = 12 — 15.5. The Hoq et al. and Shirley et al. surveys find median column 
densities of clumps drawn from the ATLASGAL and BPGS Galactic plane surveys that are 
in agreement near log(med{N}) = 14.4. Miettenen’s (2014) analysis of infrared dark clouds 
(IRDCs) indicates a significantly lower median column density of log(med{N}) = 13.1. All 
of these surveys use observations of H 13 CO + to correct the optical depth of the HCO + line. 
The Barnes et al. (2011) blind mapping survey of a section of the southern Galactic plane 
find a smaller value of log(med{N}) = 12.9, but did not observe H 13 CO + to make optical 
depth corrections. We chose a column density that is intermediate between the IRDC, low- 
mass samples, and the Galactic plane surveys of log(We/) = 14.0. For comparison, Evans 
(1999) tabulate the effective density for HCO + assuming log A" = 13.5. 

The reference column densities for each molecule in Table 1 are estimated from surveys in 
dense gas tracers (Suzuki et al. 1992, Zylka et al. 1992, Mangum & Wootten 1993, Plume et 
al. 1997, Jijina et al. 1999, Savage et al. 2002, Pirogov et al. 2003, Purcell et al. 2006, Blair 
et al. 2008, Rosolowsky et al. 2008, Falgarone et al. 2008, Hily-Blant et al. 2010, Barnes et 
al. 2011, Dunham et al. 2011, Ginsburg et al. 2011, Reiter et al. 2011, Adande et al. 2012, 
Benedettini et al. 2012, Wienen et al. 2012, Hoq et al. 2013, Shirley et al. 2013, Gerner et 
al. 2014, Miettenen et al. 2014 Yuan et al. 2014). Since these surveys contain many biases 
(i.e. only study massive high-mass regions or only study nearby low-mass regions, etc.), N re f 
represents an educated guess at a typical value intermediate between high-mass clumps and 
low-mass cores. The tabulated n e f / can be easily scaled to another column density using 
Equation 13. 13 C isotopologue column densities assume no fractionation and an isotope ratio 
of [ 12 C]/[ 13 C] = 50, appropriate for a source in the Molecular Ring (Wilson & Rood 1994). 
For molecules with ortho and para states, LTE statistical equilibrium was assumed (2 : 1 
ortho to para ratio for NH 3 and 3 : 1 ortho to para ratio for H 2 CO). 
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3.2. Comparing Critical Density and Effective Excitation Density 

The panels of Figure 4 summarize the range in n^' no bg and n e // of commonly observed 
molecular tracers of dense gas for reference observed column densities and a gas kinetic 
temperature of T k = 15 K. This temperature was chosen to be intermediate between the 
typical temperature of nearby low-mass cores (10 K, Jijina et al. 1999, Rosolowsky et al. 
2008, Seo et al. 2015, submitted) and high-mass clumps in the Galactic plane (20 K, Dunham 
et al. 2011, Wienen et al. 2012). n e j / is lower than n^ ,nob9 in nearly all cases by 1 — 2 orders 
of magnitude. This is partially due to the sensitivity of n e // to line trapping in optically 
thick transitions. Transitions that are the most optically thick have the largest average 
r ^cru ' n0 b °/ n eff ratios (be 121 for HCO + transitions and 88 for HCN transitions where the 
average is calculated for 1 — 0 through 4 — 3 transitions). Less optically thick species have 
smaller ratios (i.e. average ratios of 16 for CN and 10 for N 2 H + ). The species plotted with 
the smallest difference between n^’ nob9 and n e // is the optically thin isotpologue H 13 CN 
with an average ratio of 1.7. Some heavier species whose transitions are likely close to 
optically thin listed in Table 1 also have low average ratios that are less than 3 (i.e. H 2 CO, 
HC 3 N, and CH 3 CN). These results indicate that a 1 K km/s line is typically observed for 
transitions when 3? > 1 in Equation 1 or when the radiative rates are greater than collisional 
depopulation out of the upper energy level of the transition. 

There is also typically an increase in bg and n e ff with J u for a particular molecule 

(Figure 4). This behavior can be understood from the functional dependences on A jk and 7 j k 
in Equation 4 {n^’ noh9 ~ Ajk/'jjk)- The Einstein A for electric dipole allowed transitions 
is Ajk oc nlvj k , where fi e is the molecular permanent electric dipole moment. In general, 
molecules with larger /i e will trace denser gas. For instance, /j e = 2.99 D for HCN while 
/j, e = 0.11 D for CO (Ebenstein & Muenter 1984, Goorvitch 1994), partially explaining why 
HCN is a denser gas tracer than CO (the other effect being that the optical trapping in 
12 CO lines is often very large; see §2.2). Since vj k oc J u and the Einstein A jk has a 
dependence, higher frequency transitions of a particular molecule will trace higher densities. 
This statement is also true for comparing molecules with similar permanent electric dipole 
moments but with transitions at different frequencies. As an example, the frequency of the 
NH 3 (/i e = 1.47 D; Ueda & Iwahori 1986) inversion transitions (J,K=J) are nearly a factor 
of 4 lower than for the 1 — 0 transitions of HCO + N 2 H + and HCN, thus the spontaneous 
transition rate is more than 60 times smaller, lowering the critical densities of NH 3 (1,1) 
with respect to the 1 — 0 transitions of HCO + N 2 H + and HCN. One problem with the use of 
becomes apparent for the NH 3 inversion transitions. Since NH 3 (1,1) through (4,4) are 
at nearly the same frequency (23.7 to 24.1 GHz), their n^’ nob9 are nearly identical. n e // 
is sensitive to both line trapping and E u /k of the transition resulting in an increase in n e f / 
for NH 3 (1,1) and (2,2) (N.B. n e ff for NH 3 (3,3) and (4,4) are undefined at T k = 15 K as we 
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will discuss below). 

Differences in collision rates also affect the ranking of molecules with density in Figure 
4. HCN is observed to have 1-2 orders of magnitude higher n* , r l *"’ no hg and n e ff than HCO + . 
Both HCO + and HCN transitions have similar Aj k values, but HCO + is an ion and the 
collisional cross section for collisions with H 2 is larger than for neutral species (see Flower 
1999 vs. Dumouchel et al. 2010). This reduces the critical density and effective excitation 
density of HCO + relative to HCN. 

The ordering of species with increasing density is different for n tbl ^' noba versus n e f / 
among many of the molecules plotted in Figure 4. A problem with using n ^™’ 110 bfi can be 
seen when HCO + and N 2 H + are compared. HCO + and N 2 H + are both ions with similar A ]k 
values for transitions and similar collisional cross-sections and thus have similar b °\ 

however, it is common knowledge to observers of nearby, low-mass dense cores that N 2 H + is a 
denser gas tracer than HCO + . Due to chemical effects, N 2 H + typically has lower abundance 
and therefore lower column density than HCO + resulting in a higher n e f / (see Equation 13). 

n e ff also does a better job than n t ^ ,nobg of characterizing the difference in density 
traced by isotopologues versus the main species. The ratio of Aj k for the 1 —> 0 transition 
between H 13 CN and HCN is 0.92. This accounts for the slightly lower n t ^ ,nob9 for H 13 CN 
versus HCN; however, we know that the isotopologues should trace higher density gas than 
the main species because it is more optically thin and has a lower column density (see 
Equation 13). The average ratio of n e ff for H 13 CN/HCN 1 —> 0 is 41 ( T k e [10,100] K). 
This ratio is not exactly equal to the assumed isotope ratio and column density ratio of 50 
because of differences in Aj k , collisional rates, and line trapping between the isotopologue 
and main species. 

These examples illustrates important advantages in using n e ff rather than the critical 
density; however, using n e f / has a significant disadvantage in that n e // is undefined when 
an integrated intensity of 1 K krn/s cannot be observed. This is the case for the NH 3 
(3,3) ( E u /k = 123.5 K) and higher inversion lines at low gas kinetic temperature. When 
E u /k ^ T k , the upper energy state cannot attain a high enough population to result in large 
enough radiative rates to produce a 1 K krn/s line. High J u transitions for species such as 
NH 3 , HCN, HCO + , and N 2 H + will be undefined at low gas kinetic temperatures (T k < 20 
K). 

Another problem with the standard definition of n e ff is that the more complex the 
spectrum of a molecule (e.g. a heavy asymmetric top molecule), the less likely a particular 
transition will attain enough level population to produce a 1 K krn/s emission line. For 
instance, no transitions of HNCO have intensities greater than 1 K km/s for the observed 
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Fig. 4.— Left: n^ ,nob9 calculated at 27 = 15 K for popular dense gas tracers. Each point 
corresponds to the upper state rotational quantum number J of the transition. The ordering 
of labels correspond to the ordering of the curves in increasing density for each molecule. 
Note that the curves for HCO" 1 " and N 2 H + nearly overlap as well as for CN and H 13 CN. 
Right: n e ff calculated for 27 = 15 K and at the reference column density determined from 
surveys of Galactic sources (see Table 1). Note that the curves for HCN and N 2 H + nearly 
overlap. The critical density is typically one to two orders of magnitude above n e // for these 
tracers. 
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median column density of logmed{N} = 12.6 (Gerner et al. 2014). n e f / for complex organic 
species (i.e. CH 3 CHO, CH 3 OCH 3 , etc.) will be undefined. One solution is to reduce the 1 K 
km/s criterion, Ultimately, the choice for this criterion is arbitrary, but unless all researchers 
agree to use the same criterion, n e ff could not be compared among different papers and 
for different species. The upper right panel of Figure 3 shows the non-linear sensitivity of 
the integrated intensity for different transitions of HCO + with density. As the integrated 
intensity criterion is lowered, the curves flatten with n c indicating that a smaller decrease 
in the criterion results in a larger decrease in n e ff. These problems limit the utility of n e ff 
compared to n^’ no ba which can always be calculated if collision rates are known. 


4. A Caveat About Continuum Radiation Backgrounds 


We ignored continuum background radiation fields in §2 when deriving the expressions 
for ri^' no bg and n^ k ' no 6fl , while only the CMB was considered as a background in the 
RADEX models in §3. In this section we consider when contributions from the CMB and 
dust continuum background are important. 

The ratio of the rate of stimulated emission from j —> k due to the CMB to spontaneous 
emission is given by 


Bjk'U'Vjk (Tcmb) 


A 


jk 


ghvj fc / kT crn b 


^ph (T C mbi G k ) 


(14) 


This ratio is 2 for the NH 3 (J,K=J) inversion transitions and is a modest 0.26 at 90 GHz 
(3.3 mm), near the frequency of the fundamental transitions of commonly observed dense 
gas tracers (e.g. HCO + , HCN, HNC, N 2 H + ). This ratio drops below 0.1 for transitions 
with Ujk > 140 GHz (A jk < 2.1 mm). Table 1 lists n p h(T C mb, Vjk) for each transition. For 
low frequency (long wavelength) transitions, stimulated emission and absorption due to the 
CMB should not be ignored. 


Interesting examples of the interaction of low frequency transitions with the CMB are 
the K-doublet transitions of ortho-H 2 CO (l^o — lpi, 2^1 — 2 1j2 , etc.). The 4.6 GHz (6.5 cm) 
li,o — Iqi transition has a stimulated emission rate to spontaneous emission rate ratio of 11.3. 
Since the statistical weights of the two level are identical, the absorption rate from the lower 
level is identical to the stimulated emission rate from the upper level. This line is typically 
observed in absorption for densities n < 10 5 enr 3 (Mangum et al. 2008). The absorption is 
created by the lowest energy level (lgi) having a collisionally-induced overpopulation that 
results in absorption against the CMB ( T ex < T cmb ) over a wide range of physical conditions 
(Evans 1975, Garrison 1975). The calculation of n tb r ™ ,nobg gives a ridiculously low value of 
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13 cm -3 at Tk = 10 K. Clearly in cases such as these, the simple definition of critical density 
fails and more sophisticated radiative transfer and statistical equilibrium calculations that 
include the CMB background are required. 

In addition to the CMB, dust continuum backgrounds from the Galactic interstellar 
radiation field and from emission local to the cloud can be important in modifying the level 
populations. The main effect of dust continuum emission on our definition of critical density 
is to introduce a new stimulated rate out of level j to higher energy levels, i. The energy 
density is now the energy density for optically thin dust emission 


_ dust (rj~\ \ dust (rp \ 

U'vjn \-L dust) ' i/A • ^ dust) ? 


(15) 


where r^ ust is the dust optical depth at the frequency of the i —> j transition and T dust is 
the isothermal dust temperature (see Equation 4 of Shirley et al. 2011 for a derivation of 
the equivalent isothermal dust temperature at a given frequency in sources with density and 
temperature gradients). We only need to consider transitions between i —> j that are allowed 
by electric dipole selection rules because the Einstein Bji for absorption of j —> i is directly 
proportional to the Einstein A tJ for spontaneous emission of i —> j (see Equation 2). For 
many simple molecules (HCO + , N 2 H + , HCN, CS, NH 3 , etc.) there is only one electric dipole 
permitted transition from j into the upper state i that needs to be considered. The ratio of 
the rate at which this absorption from j —> i occurs to the spontaneous rate for the original 
j —y k transition is 


r> q . dust 

Ajk 


9i a. 


V 


T, 


dust 


gj A jk exp (huij/kT dust ) - 1 

” a ^H^H 2 ^Vjk^ph (Tdust 5 l^ij ) 

9j Ajk 


(16) 


where fi = 2.8 is the mean mass per H 2 molecule (see Kauffmann et al. 2008), Nu 2 is the H 2 
column density (cm -2 ), and n Ujk is dust opacity (cm 2 /gram of gas) at the frequency of the 
j —> % absorption. The ratio of rates in Equation 16 is linear proportional to the H 2 column 
density and increases non-linearly with dust temperature. 


For the following examples, we calculate the ratio of rates in Equation 16 assuming a H 2 
column density of 10 22 cm -3 (corresponding to Ay ~ 10 mag), a gas mass to dust mass ratio 
of 100 : 1, and Ossenkopf & Henning opacities for coagulated grains with thin ice mantles 
(OH5, Ossenkopf & Henning 1994) that are appropriate for dense cores. For the NH 3 (1,1) 
transition, the upper level is radiatively coupled to the J = 2,K = 1 level via a 1.168 THz 
(267 /jrn) far-infrared transition (see Figure 1). If NH 3 is bathed in a dust continuum of 
10 K, the stimulated rate of absorption out of the upper inversion level is four times the 
spontaneous rate of the (1,1) transition. This high ratio is partially driven by the ratio of 
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spontaneous rates between the terahertz transition and the centimeter inversion transition 
which is Aij/Ajk = 7.2 x 10 5 . The net effect is that n e ff becomes larger as the energy 
density from dust continuum increases. When the upper level of a low frequency transition 
is radiativcly coupled to far-infrared transitions, the absorption from dust continuum is an 
important rate that should not be ignored. This is the case for the centimeter transitions of 
molecules such as NH 3 , OH, and CH. 

The effects of dust continuum radiation are significantly less for the lowest energy tran¬ 
sitions of simple linear molecules. For the 267.5 GHz (1.1 mm) HCO + 3 — 2 transition at 
10 K, the ratio of rates is 6 x 10~ 4 . At 50 K, the ratio of rates rises to a meager 7 x 10~ 3 . 
These very low ratios indicate that the absorption due to dust continuum emission is not 
important for the level populations associated with millimeter transitions of linear molecules 
such as HCO+, HCN, CS, and N 2 H+. 

We have only considered the rotational levels of the molecule, but there can also be 
coupling between the radiation field and vibrational levels or electronic levels of the molecule 
that can effect the level populations of the ground vibrational state rotational levels. A mid- 
infrared continuum can pump the molecule from the ground vibrational state into an excited 
vibrational state which then radiatively decays via ro-vibrational transitions. For example, 
the HCN and HNC molecules have bending modes (iq, i/ 2 ,z /3 = 0,1,0) that can be excited 
by mid-infrared pumping at 14 /irn and 21 /iin respectively. In environments where such 
pumping may be important (i.e. the inner envelope of high-mass protostars or near an 
AGN), the level populations can be affected by mid-infrared pumping. In these cases, more 
sophisticated radiative transfer models are needed. 


5. Summary 

In this tutorial, n tb ™’ nobg and n e ff are defined and their properties analyzed. The values 
for ri bl ^'' no hfJ and n e ff are calculated for multiple transitions of 12 commonly observed dense 
molecular gas species in Table 1. n e ff inversely depends on the choice of column density 
and increases non-linearly as gas kinetic temperatures are lowered (77 < 20 K). We have 
estimated the reference column density for difference species listed in Table 1 from surveys 
of cores and clumps in the Milky Way. The n e ff quoted can be easily scaled to another 
column density using Equation 13. 

Since the definition of n e // is arbitrary (the 1 K km/s criteria) and does not depend 
on a specific ratio of radiative rates to collisional depopulation rates (c.f. 37 = 1 for the 
definition of n crit ), it is seductive to use n tbl £’ nobg to characterize the density at which a 
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transition is excited as it can be analytically calculated; however, n£.™’ no bg does not account 
for radiative trapping which can easily lower the effective critical density by more than an 
order of magnitude for commonly observed dense gas tracers such as HCO + 1 — 0, HCN 
1 — 0, N 2 H+ 1 — 0, etc. The combination of 9ft > 1 and radiative trapping that generate a 1 K 
km/s line lower n e f / by ~ 1 — 2 orders of magnitude relative to n^' no bg for many commonly 
observed millimeter and submillimeter lines. Furthermore, n e ff can account for abundance 
or column density differences to better differentiate the excitation density of different species 
than n^ ,nobg (i.e. HCO + vs. N 2 H + or isotpologues vs. the main species). Overall, the 
effective excitation density gives a better estimate of the range of densities at which a modest 
(1 K km/s) molecular line can be observed; however, use of n e ff has the disadvantages that 
an appropriate column density must be determined for the observations of interest and that 
it is not defined for all transitions (in particular when E u /k 3> T*,) and for all species (e.g. 
HNCO, CH 3 CHO, etc.). 

n cr1" n0 bg an d n eff provide only rough estimates of the densities traced and should not be 
over-interpreted. More sophisticated tools such as the Contribution Function (Tafalla et al. 
2006) may be used to determine the various contributions to the observed line profile along 
the line of sight (see Pavlyuchenkov et al. 2008 for a detailed analysis of these techniques). 
Ultimately, if one wants to understand their observed spectra, radiative transfer modeling 
with publicly available codes are a fast and efficient way to determine the physical properties 
of a region. 
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Table 1. and n e // for Dense Gas Tracers with Collisions with H 2 


Molecule 

3 -> k 

u jk 

(GHz) 

Bj/k 

(K) 

Ajk 

(s- 1 ) 

n ph (.T C mb ) 


thin.no bq 
n crit 

( T k ) cm 

-3 

log N re f 
(cm -2 ) 


n eff 

(T k ,N ref )* 

cm -3 


Ref. 

10K 

20K 

50K 

100K 

10K 

15K 

20K 

50K 

100K 

HCO+ 

1-0 

89.189 

4.28 

4.3E-5 

0.264 

6.8E+4 

4.5E+4 

2.9E+4 

2.3E+4 

14.0 

9.5E+2 

6.4E+2 

5.3E+2 

3.3E+2 

2.6E+2 

1 


2-1 

178.375 

12.84 

4.1E-4 

0.046 

5.6E+5 

4.2E+5 

2.8E+5 

2.2E+5 


4.1E+3 

2.3E+3 

1.7E+3 

8.9E+2 

6.5E+2 



3-2 

267.558 

25.68 

1.5E-3 

0.009 

1.6E+6 

1.4E+6 

1.0E+6 

8.1E+5 


2.4E+4 

1.1E+4 

6.8E+3 

2.6E+3 

1.7E+3 



4-3 

356.734 

42.80 

3.6E-3 

0.002 

3.6E+6 

3.2E+6 

2.5E+6 

2.0E+6 


1.9E+5 

5.4E+4 

2.8E+4 

7.2E+3 

4.0E+3 


h 13 co+ 

1-0 

86.754 

4.16 

3.9E-5 

0.279 

6.2E+4 

4.1E+4 

2.7E+4 

2.0E+4 

12.3 

3.9E+4 

2.7E+4 

2.2E+4 

1.4E+4 

1.1E+4 

1 


2-1 

173.507 

12.49 

3.7E-4 

0.050 

5.1E+5 

3.8E+5 

2.6E+0 

2.0E+5 


1.7E+5 

9.6E+4 

7.1E+4 

3.7E+4 

2.6E+4 



3-2 

260.255 

24.98 

1.3E-3 

0.011 

1.5E+6 

1.3E+6 

9.5E+5 

7.3E+5 


1.3E+6 

4.3E+5 

2.7E+5 

9.6E+4 

6.2E+4 



4-3 

346.998 

41.63 

3.3E-3 

0.002 

3.4E+6 

2.9E+6 

2.3E+6 

1.8E+6 



2.7E+6 

1.1E+6 

2.7E+5 

1.6E+5 


n 2 h+ 

1-0 

93.174 

4.47 

3.6E-5 

0.242 

6.1E+4 

4.1E+4 

2.6E+4 

2.0E+4 

13.0 

1.0E+4 

6.7E+3 

5.5E+3 

3.4E+3 

2.6E+3 

1,2 


2-1 

186.345 

13.41 

3.5E-4 

0.040 

5.0E+5 

3.7E+5 

2.6E+5 

1.9E+5 


4.6E+4 

2.6E+4 

1.9E+4 

9.7E+3 

6.9E+3 



3-2 

279.512 

26.83 

1.3E-3 

0.007 

1.4E+6 

1.2E+6 

9.2E+5 

7.1E+5 


2.6E+5 

1.1E+5 

6.8E+4 

2.5E+4 

1.6E+4 



4-3 

372.673 

44.71 

3.1E-3 

0.001 

3.2E+6 

2.8E+6 

2.2E+6 

1.7E+6 


3.0E+6 

5.4E+5 

2.7E+5 

6.3E+4 

3.4E+4 


HCN 

1-0 

88.632 

4.25 

2.4E-5 

0.268 

4.7E+5 

3.0E+5 

1.7E+5 

1.1E+5 

14.0 

8.4E+3 

5.6E+3 

4.5E+3 

2.6E+3 

1.7E+3 

3 


2-1 

177.261 

12.76 

2.3E-4 

0.047 

4.1E+6 

2.8E+6 

1.6E+6 

1.1E+6 


4.1E+4 

2.2E+4 

1.6E+4 

7.3E+3 

4.5E+3 



3-2 

265.886 

25.52 

8.4E-4 

0.010 

1.4E+7 

1.0E+7 

5.7E+6 

3.8E+6 


3.0E+5 

1.2E+5 

7.3E+4 

2.5E+4 

1.3E+4 



4-3 

354.505 

42.53 

2.1E-3 

0.002 

3.0E+7 

2.3E+7 

1.4E+7 

9.1E+6 


2.1E+6 

5.8E+5 

3.2E+5 

8.1E+4 

3.7E+4 


h 13 cn 

1-0 

86.340 

4.14 

2.2E-5 

0.282 

5.3E+5 

2.5E+5 

1.3E+5 

9.7E+4 

12.3 

3.5E+5 

2.2E+5 

1.6E+5 

8.7E+4 

6.5E+4 

4 


2-1 

172.678 

12.43 

2.1E-4 

0.051 

3.4E+6 

2.2E+6 

1.2E+6 

9.1E+5 


1.9E+6 

8.6E+5 

5.8E+5 

1.8E+5 

1.7E+5 



3-2 

259.012 

24.86 

7.7E-4 

0.011 

8.8E+6 

6.6E+6 

4.1E+6 

3.3E+6 


2.5E+7 

3.8E+6 

2.0E+6 

6.3E+5 

4.2E+5 

1 


4-3 

345.340 

41.43 

1.9E-3 

0.002 

1.9E+7 

1.5E+7 

9.7E+6 

7.7E+6 



5.6E+7 

1.1E+7 

1.9E+6 

1.1E+6 

to 

HNC 

1-0 

90.664 

4.35 

2.7E-5 

0.256 

1.4E+5 

1.1E+5 

8.4E+4 

7.0E+4 

14.0 

3.7E+3 

2.7E+3 

2.3E+3 

1.7E+3 

1.3E+3 

3 ^ 


2-1 

181.325 

13.05 

2.6E-4 

0.043 

1.3E+6 

1.0E+6 

8.2E+5 

6.8E+5 


2.0E+4 

1.2E+4 

9.3E+3 

5.5E+3 

3.9E+3 

1 


3-2 

271.981 

26.11 

9.3E-4 

0.009 

5.1E+6 

4.0E+6 

3.1E+6 

2.5E+6 


1.6E+5 

7.4E+4 

4.9E+4 

2.1E+4 

1.2E+4 



4-3 

362.630 

43.51 

2.3E-3 

0.002 

1.3E+7 

1.0E+7 

7.8E+6 

6.2E+6 


1.4E+6 

4.2E+5 

2.4E+5 

6.8E+4 

3.4E+4 


CN 

1 3/2 — °l/2 

113.495 

5.45 

1.2E-5 

0.158 

4.1E+5 

2.4E+5 

1.1E+5 

6.4E+4 

14.0 

3.8E+4 

2.3E+4 

1.7E+4 

7.7E+3 

4.6E+3 

5 


2 5/2 — 1 3/2 

226.876 

16.34 

1.1E-4 

0.019 

3.2E+6 

2.2E+6 

1.1E+6 

6.4E+5 


2.3E+5 

1.1E+5 

7.3E+4 

2.8E+4 

1.5E+4 



37/2 — 2 5/2 

340.249 

32.66 

4.1E-4 

0.003 

8.3E+6 

6.5E+6 

3.5E+6 

2.1E+6 


1.8E+6 

5.6E+5 

3.6E+5 

1.1E+5 

5.1E+4 


CS 

1-0 

48.991 

2.40 

1.7E-6 

0.734 

1.5E+4 

1.1E+4 

7.8E+3 

5.7E+3 

13.5 

8.7E+3 

6.4E+3 

5.4E+3 

3.4E+3 

2.4E+3 

6 


2-1 

97.981 

7.10 

1.7E-5 

0.218 

1.3E+5 

1.0E+5 

7.4E+4 

5.5E+4 


2.2E+4 

1.5E+4 

1.2E+4 

6.9E+3 

4.7E+3 



3-2 

146.969 

14.10 

6.1E-5 

0.082 

4.4E+5 

3.6E+5 

2.6E+5 

2.0E+5 


7.4E+4 

4.4E+4 

3.3E+4 

1.6E+4 

1.0E+4 



4-3 

195.954 

23.50 

1.5E-4 

0.033 

1.0E+6 

8.6E+5 

6.4E+5 

4.8E+5 


2.8E+5 

1.4E+5 

9.2E+4 

3.6E+4 

2.0E+4 



5-4 

244.936 

35.30 

3.0E-4 

0.014 

2.0E+6 

1.7E+6 

1.3E+6 

9.5E+5 


1.4E+6 

4.2E+5 

2.5E+5 

7.6E+4 

3.8E+4 



6-5 

293.912 

49.40 

5.2E-4 

0.006 

3.4E+6 

2.9E+6 

2.2E+6 

1.6E+6 



1.5E+6 

7.1E+5 

1.6E+5 

7.2E+4 



7-6 

342.883 

65.80 

8.4E-4 

0.002 

5.0E+6 

4.4E+6 

3.4E+6 

2.6E+6 



8.8E+6 

2.1E+6 

3.4E+5 

1.3E+5 


p-NH 3 c 

(1,1) 

23.694 

1.14 

1.7E-7 

1.940 

2.0E+3 

1.8E+3 

1.2E+3 

8.7E+2 

14.3 

9.3E+2 

7.9E+2 

7.6E+2 

7.4E+2 

7.3E+2 

7 


(2,2) 

23.723 

42.32 

2.3E-7 

1.937 

1.6E+3 

1.5E+3 

1.4E+3 

1.1E+3 



1.6E+4 

3.9E+3 

1.3E+3 

9.4E+2 



2+ _ i — 

Z 1 1 l_ 

1168.452 

57.21 

1.2E-2 

0.000 

1.2E+8 

1.1E+8 

9.0E+7 

6.7E+7 



4.8E+7 

8.6E+6 

1.4E+6 

8.2E+5 



2i -1+ 

1215.245 

58.32 

1.4E-2 

0.000 

1.4E+8 

1.2E+8 

1.0E+8 

7.5E+7 



7.2E+7 

1.0E+7 

1.6E+6 

8.9E+5 


o-NH 3 c 

(3,3) 

23.870 

123.54 

2.6E-7 

1.923 

1.5E+3 

1.5E+3 

1.4E+3 

1.2E+3 

14.6 




4.0E+2 

1.7E+2 

7 


0 

1 

0 

0 

572.498 

27.48 

1.6E-3 

0.000 

4.0E+7 

3.1E+7 

1.7E+7 

9.9E+6 


6.2E+5 

1.5E+5 

8.4E+4 

3.1E+4 

2.5E+4 



to 

0 

1 

0 

1214.859 

85.78 

1.8E-2 

0.000 

2.1E+8 

1.9E+8 

1.6E+8 

1.1E+8 



4.2E+7 

4.4E+6 

5.2E+5 

3.0E+5 


p-H 2 CO 

lo,l — °0,0 

72.838 

3.50 

8.2E-6 

0.386 

4.5E+4 

2.8E+4 

1.8E+4 

1.4E+4 

12.7 

5.0E+4 

3.2E+4 

2.6E+4 

1.9E+4 

1.8E+4 

8 


2(3,2 — 1(3,1 

145.603 

10.50 

7.8E-5 

0.084 

3.4E+5 

2.5E+5 

1.7E+5 

1.3E+5 


1.5E+5 

8.2E+4 

6.3E+4 

4.3E+4 

4.0E+4 



3(3,3 — 2(3,2 

218.222 

21.00 

2.8E-4 

0.022 

9.7E+5 

7.8E+5 

5.7E+5 

4.7E+5 


7.7E+5 

2.9E+5 

2.0E+5 

1.2E+5 

1.1E+5 



4(3,4 — 3o,3 

290.623 

34.90 

6.9E-4 

0.006 

1.9E+6 

1.7E+6 

1.3E+6 

1.1E+6 



1.5E+6 

7.3E+5 

3.1E+5 

2.6E+5 



5(3,5 — 4(3,4 

362.736 

52.30 

1.4E-3 

0.002 

3.6E+6 

3.1E+6 

2.6E+6 

2.2E+6 




3.8E+6 

8.0E+5 

5.9E+5 


o-H 2 CO 

ll,0 — li,i 

4.830 

15.40 

3.6E-9 

11.310 

1.3E+1 

1.0E+1 

7.6E+0 

6.1E+0 

13.0 






8 


2 l,l — 2 i )2 

14.488 

22.60 

3.2E-8 

3.456 

1.1E+2 

8.7E+1 

6.5E+1 

5.3E+1 









2l,2 — ll,l 

140.840 

21.90 

5.3E-5 

0.092 

1.5E+5 

1.3E+5 

1.1E+5 

9.0E+4 


9.5E+4 

5.9E+4 

4.6E+4 

2.6E+4 

2.0E+4 



2 i,i — li,o 

150.498 

22.60 

6.5E-5 

0.077 

2.2E+5 

1.8E+5 

1.3E+5 

1.1E+5 


1.7E+5 

1.1E+5 

8.6E+4 

4.9E+4 

3.5E+4 




Table 1—Continued 


Molecule 

i -> k 

u jk 

(GHz) 

Ej/ k 
(K) 

Ajk 

(s- 1 ) 

n ph( T cmb) 


thin,no bg 
n crit 

(T fe ) cm 

-3 

log N re f 
(cm" 2 ) 


n eff 

(■ T k ,N ref r 

cm -3 


Ref. 

10K 

20K 

50K 

100K 

10K 

15K 

20K 

50K 

100K 


3l,3 — 2 1,2 

211.211 

32.10 

2.3E-4 

0.025 

5.6E+5 

5.1E+5 

4.3E+5 

3.8E+5 


4.6E+5 

2.0E+5 

1.4E+5 

6.3E+4 

4.9E+4 



3i,2 — 2 1,1 

225.698 

33.40 

2.8E-4 

0.019 

8.4E+5 

7.0E+5 

5.5E+5 

4.6E+5 


8.9E+5 

3.8E+5 

2.6E+5 

1.3E+5 

9.0E+4 



4 1,4 — 3i,3 

281.527 

45.60 

5.9E-4 

0.007 

1.4E+6 

1.3E+6 

1.1E+6 

9.8E+5 


2.9E+7 

9.9E+5 

5.1E+5 

1.7E+5 

1.2E+5 



4 1,3 — 3l,2 

300.837 

47.90 

7.2E-4 

0.005 

2.0E+6 

1.7E+6 

1.4E+6 

1.2E+6 



1.9E+6 

9.6E+5 

3.3E+5 

2.2E+5 



5 1,5 — ^1,4 

351.769 

62.50 

1.2E-3 

0.002 

2.8E+6 

2.6E+6 

2.2E+6 

2.0E+6 



1.7E+7 

2.3E+6 

4.3E+5 

2.6E+5 


HC 3 N d 

2-1 

18.196 

1.31 

3.9E-7 

2.659 

9.7E+2 

7.5E+2 

5.3E+2 

4.1E+2 

13.0 

5.9E+3 

3.1E+3 

2.3E+3 

1.0E+3 

8.4E+2 

9 


3-2 

27.294 

2.62 

1.4E-6 

1.628 

3.5E+3 

2.6E+3 

1.8E+3 

1.6E+3 


4.7E+3 

3.0E+3 

2.3E+3 

1.1E+3 

8.7E+2 



4-3 

36.392 

4.37 

3.5E-6 

1.119 

9.1E+3 

6.8E+3 

4.7E+3 

3.7E+3 


6.2E+3 

3.8E+3 

2.8E+3 

1.5E+3 

1.1E+3 



5-4 

45.490 

6.55 

6.9E-6 

0.819 

1.8E+4 

1.3E+4 

9.6E+3 

7.2E+3 


9.6E+3 

5.6E+3 

4.1E+3 

2.1E+3 

1.5E+3 



8-7 

72.784 

15.72 

2.9E-5 

0.387 

7.7E+4 

5.9E+4 

4.1E+4 

3.1E+4 


5.3E+4 

2.5E+4 

1.6E+4 

6.4E+3 

4.0E+3 



9-8 

81.881 

19.65 

4.2E-5 

0.312 

1.1E+5 

8.3E+4 

6.2E+4 

4.6E+4 


1.1E+5 

4.2E+4 

2.6E+4 

9.4E+3 

5.7E+3 



10-9 

90.979 

24.01 

5.8E-5 

0.254 

1.6E+5 

1.2E+5 

8.2E+4 

6.3E+4 


4.3E+5 

7.2E+4 

4.3E+4 

1.4E+4 

8.1E+3 



11 - 10 

100.076 

28.82 

7.8E-5 

0.209 

2.0E+5 

1.6E+5 

1.1E+5 

9.2E+4 



1.6E+5 

7.0E+4 

2.0E+4 

1.1E+4 



12 - 11 

109.173 

34.06 

1.0E-4 

0.173 

2.8E+5 

2.1E+5 

1.5E+5 

1.1E+5 



2.3E+5 

1.1E+5 

2.8E+4 

1.6E+4 



14 - 13 

127.367 

45.85 

1.6E-4 

0.120 

3.4E+5 

3.1E+5 

2.5E+5 

1.9E+5 




3.8E+5 

5.7E+4 

3.0E+4 



15 - 14 

136.464 

52.40 

2.0E-4 

0.100 

5.5E+5 

4.3E+5 

3.1E+5 

2.4E+5 




8.2E+5 

8.1E+4 

4.2E+4 


CH 3 CN e 

4 0 - 3 0 

73.590 

8.83 

3.2E-5 

0.379 

1.2E+5 

8.6E+4 

5.4E+4 

4.1E+4 

13.0 

5.2E+4 

4.0E+4 

3.7E+4 

2.5E+4 

2.5E+4 

10 


5o - 4 q 

91.987 

13.24 

6.3E-5 

0.249 

2.5E+5 

1.7E+5 

1.1E+5 

8.3E+4 


1.2E+5 

8.4E+4 

7.4E+4 

4.3E+4 

4.0E+4 



6 0 - 5 0 

110.384 

18.54 

1.1E-4 

0.169 

4.5E+5 

3.1E+5 

1.9E+5 

1.5E+5 


3.2E+5 

1.8E+5 

1.5E+5 

7.6E+4 

6.5E+4 



7q - 6 0 

128.779 

24.72 

1.8E-4 

0.117 

7.3E+5 

5.1E+5 

3.1E+5 

2.4E+5 


1.0E+6 

4.2E+5 

3.1E+5 

1.4E+5 

1.1E+5 



00 

0 

1 

0 

147.175 

31.79 

2.7E-4 

0.082 

1.1E+6 

7.7E+5 

4.8E+5 

3.6E+5 


8 .8E+6 

1.0E+6 

6.5E+5 

2.3E+5 

1.7E+5 



9q — 8q 

165.569 

39.73 

3.8E-4 

0.058 

1.6E+6 

1.1E+6 

6.9E+5 

5.2E+5 



3.1E+6 

1.4E+6 

4.0E+5 

2.8E+5 



a The effective excitation density may be converted to another column density using n e ff(Tp c , N ) 
to convert ra e y y to other Tp ,. 


n eff( T k’ N ref) for n eff < 5 X 10® cm 3 . 


Logarithmic interpolation should be used 


^The collisional rates for collisions with H 2 in cm -3 s — ^ were obtained from the Leiden LAMDA database (August 2014). The original references for the collision rates are: 1. Flower et al. 
1999, 2. Daniel et al. 2005, 3. Dumouchel et al. 2010, 4. Green 8z Thaddeus 1974, 5. Lique et al. 2011, 6. Lique et al. 2006, 7. Danby et al. 1988, 8. Wiesenfeld 8z Faure 2013, 9. Green 8z 
Chapman 1978, 10. Green 1986. Rates calculated for collisions with He are scaled in the LAMBDA database by the reduced mass (1.36) to convert to approximate rates for collisions with H 2 - 
Collision rates with hyperfine transitions collapsed were used for the calculations in this table. 

Cn crit l ' Tl ° was logarithmically extrapolated to Tp~ = 10K since the lowest Tp. calculated in the collision rate file was 15K. 

C * n crii l,n ° was logarithmically extrapolated to Tp~ = 100K since the highest Tp~ calculated in the collision rate file was 80K. 


was logarithmically extrapolated to = 10K since the lowest Tp, calculated in the collision rate file was 20K. 


e thzn.no bq 
erzt 



